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Analysis of PML Method for Stochastic Convected Helmholtz 
Equation 

Sang-Hyeon Park and Imbo Sim 


Abstract. We propose and analyze the perfectly matched layer (PML) method for the time- 
harmonic acoustic waves driven by the white noise source in the presence of the uniform flow. 
A PML is an artificial absorbing layer commonly used to truncate computational regions to solve 
problems in unbounded domains. We study a modification of PML method based on Becache et. 
al. [!]• A truncated domain problem for stochastic convected Helmholtz equation in the infinite 
duct is constructed by applying PMLs. Our PML method omits the instability of inverse upstream 
modes in the PML. Moreover, a suitable jump condition on boundaries between computational 
domain and PMLs is not required. We analyze the stochastic error generated by truncations of the 
domain. Thus the convergence analysis of the solution is provided in the sense of mean-square. 
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1. Introduction 

We study the stochastic convected Helmholtz equation with the white noise source in the infinite 
duct. Let the duct be P = {a; = (a;i,a; 2 )| G R,a ;2 G (0,d)} C where d is a positive constant. 
The acoustic pressure held p(x) in the presence of a uniform how satishes the convected Helmholtz 
equation in fl: 

+ 5 + +^"p=f+ ( 1 - 1 ) 

where the Mach number M = ^, the wave number k = ■^, and cq sound velocity in the huid. In 
addition, the pressure satishes the homogeneous Neumann boundary condition on the wall, 

2^ = 0 on dP. (1.2) 

0X2 

We assume that 0 < M < 1. The right-hand side of dni) consists of deterministic term / and 
stochastic term W{x). The stochastic term W{x) is the formal derivative of the white noise in space 
ED. We suppose that the right-hand side of dm has a compact support P/ in P i.e. P/ = supp(/) 
and Py C P. Here, / G L^(P) is assumed. The function lQj,(a;) is an indicator function. 

Convected Helmholtz equation describes the time-harmonic acoustic waves generated by moving 
media. In the presence of a uniform how, time-harmonic linearized Euler equations reduce to the 
convected Helmholtz equation for the pressure. The presence of a mean how causes the mathematical 
treatment of the problem much more difficult, since phase and group velocities of the acoustic wave 
could have opposite signs [iiiiin]. 
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The main goal of this paper is to construct a modified perfectly matched layer(PML) method 
for the problem (11.11) and to analyze a suitable error estimation in the sense of i?[|| • ||| 2 ]- 

Perfectly matched layer was introduced by Berenger in [7] for the time-dependent Maxwell equa¬ 
tions. To design an efficient absorbing boundary condition for the time harmonic acoustic waves in 
unbounded domains, the PML method is usually used. In the presence of a flow, compared to the 
classical waves, PMLs could cause instabilities of the solution. The presence of instabilities have been 
studied in [5] by using group and phase velocities. In [4], Becache et. al. studied a PML formulation 
for the convected Helmholtz equation in a duct to avoid instabilities of the solution. They used a sub¬ 
stitution —>■ a-^ -I- iA. In [4], A is chosen as a constant in the PMLs. The complex-valued 

function a satisfies Re(a) > 0 and Im(a) < 0 in PMLs. This approach provides a successful analysis 
of error estimation under a suitable jump condition between the computational domain and PMLs. 

According to mm US], regularity estimates for SPDEs are usually weak which leads to low 
order error estimates. To overcome this complication, authors in [a nm HU na US] considered the 
discretized white noise instead of W, since it is more regular. Especially, Cao et. al. m established an 
error estimation of the PML method for Helmholtz equation by applying the discretized white noise. 
We also utilize the discretized white noise to obtain the desired regularity for the error estimation 
as in m- Replacing the white noise of the forcing term with a discretized version, an approximated 
problem is built as follows. 

^ = 0 Ian/ ' 

0X2 

where (x) is the discretized white noise. Since the problem (11.31) is more regular than the problem 
(HU, the convergence analysis can be carried out for the PML method as in m 

In this setting (II.3L we study a modified PML method instead of that in |4]. The modified PML 
formulation considered in this paper is described in the following. The left-hand side of the problem 
HU is reformulated by 




Mki 2 h h 


Then a substitution 

d Mki / 9 Mki 

is applied to construct a PML problem [20] . Here, the complex function a{xi) is given by 


a{xi) 


—ioj 

—iu! -I- < 7 ( 3 : 1 ) ’ 


where w > 0 and < 7 ( 3 : 1 ) is a positive real-valued function in C^(]R) such that < 7 ( 3 : 1 ) = 0 in the compu¬ 
tational domain. Therefore, is continuously connected on the boundary between the computational 
domain and PMLs. In this reason, this PML does not require a jump condition. It also follows from 
the modal expansion approach in |4] that a PML-truncation error of this model. 

The outline of this paper is as follows. In Section 2, we discuss the discretized white noise. The 
error estimation E[||p — and convergence analysis are provided. In Section 3, we discuss about 

a modified PML method and the convergence of the problem as the length of layer goes to infinity. In 
Section 3.1, problems are restricted in a bounded domain fit, in order to analyze an error generated 
by PML . In section 3.2, we derive || • ||^fi-error of the PML model. Combining two error estimations 
in section 2 and section 3.2, we finally establish the total error analysis in the sense of E[|| • \\‘j^ 2 ]- 
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2. Approximation driven by the Discretized White Noise 

In this section, we derive the approximated problem m from (HI). This approximation provides 
more regular solution than the solution of dnj. Let T be a trianglation Uti C tl/- Each element 
Ki has at most one curved edge aligned with H. Let hi = diam Ki, h = max^ hi and Vi be the radius 
of the largest ball inscribed Ki. We assume that T is quasi-uniform. The random variable ^i is define 
by 

l<*<iV. (2.1) 

Here, \Ki \ denotes the area of Ki. According to [21], the family of random variables is independent 

identically distributed (i.i.d.) normal random variables with mean 0 and variance 1. The piecewise 
constant approximation W^{x) G is defined by 


N 








( 2 . 2 ) 


which has following properties: 

E[W\x)W’^{y)] = -^, t 0 Tx,yeK,, 

|Ai| 

E[W^{x)W^(y)] = 0, for x € Ki and y G Kj, i ^ j. 

Therefore, the problem (HD can be reformulated by replacing W with W^. Consider the following 
approximated problem 

(1 _ + + 2ikM^ + tV = / + i. n, 

OX^ 0X2 (JX\ , . 

^ = 0 on dn. 

0X2 

The problem (12.31) is more regular than the problem (II.ip . since W^{x) G L^(H). Moreover, a weak 
form of the problem (12.31) is 


an{p\q) = -[{f + W^)qdx, yqeH\n), 
Jn 

where the sesquilinear form af 2 (-, •) is 


anijp.q) 


{M Vp) • Vq dx 


[ k^pqdx 

Jn 



Mp 

0 


q dx. 


(2.4) 


(2.5) 


where M = 


1 - m2 0 \ 
0 1 )■ 


Now, we prove that the solution p^ of ()2.3p converges to p in a given bounded domain Hh G H. 
Using the stochastic integration formula in [^ and [^, the mean-square error A-db — is 

derived in the following. Let G(x, y) be the Green function of the convected Helmholtz equation (12.3p . 
We denote H^\x) is a Hankel function of the first kind of order 0. By the method of images([T9j), 
for a fixed p G H, the Green function G(x, y) is written by 


G(x, y) =4>(x, y) -f 4>(x, -y) 

00 

+ XI ^n-y) + ^(2;, el + y) + $(x, e~-y) + $(x, e~ + p)), 

n—1 


( 2 . 6 ) 







4 


Sang-Hyeon Park and Imbo Sim 


where 

^ y/xj + (1 - M^)xl 

1 - M2 

et = (0,±2dn). 

Accrording to [9] and [21] , p and p^ can be written by 

p{x) = f G{x,y)f{y) dy + f G{x,y)lQ^{y) dW{y), 

Jn Jn 

p’^{x)= [ G{x,y)f{y) dy + [ G{x,y)lnj:{y) dW'^{y). 

Jn Jn 

To analyze the error with respect to p and p^, we start with following two lemmas. 


(2.7) 


( 2 . 8 ) 


Lemma 2.1. There exists a Lipschitz continuous function V(x, y) with respect to x and y such that 
$(a:, y) is decomposed by 

Proof. Let us define a linear map T : E,^ —^ by T{x) = Let x = T{x) and 

y = T{y). Then it follows from (12.71) that 

p{x - y) = \x - y\. (2.10) 


It is well-known that there exists a Lipschitz continuous function V(x, y) with respect to x and y such 
that 

iH^{k\x-y\) = ^\nj^^-^^+V{x,y) (2.11) 

(refer to [11]). Since the map T is linear, V{x,y) := V(T{x),T{x))e~^ K^i-vi) jg g^jj Lipschitz 
continuous with respect to x and y. Consequently, the proof is completed from (12.711 and (12.111) . □ 


From dill) and Lemma lO the next Lemma is proved. 

Lemma 2.2. Let Lib be an open bounded domain in R^. IFe assume that Llf G Lib G LI. Suppose 
y, z G Ki and h is sufficiently small. Then we have 


|G(x, y) - G{x, z)fdx < C\y - z|' 


( 2 . 12 ) 


where 0 < e < 1 and G is a positive constant independent of e. 

Proof. Refer to the proof in the appendix. □ 


Applying Lemma 12.21 we are now in position to obtain the error estimation 

Theorem 2.3. Let p and p^ are the solution of fny, El, and El, respectively. Let Lib be an open 
bounded domain in R^. We assume that Llf G Lib G LI. Suppose y,z G Ki and h is sufficiently small. 
Then we obtain 

i^[lb-/lli.(n,)] < (2.13) 

where 0 < e < 1 and G is a positive constant independent of e and h. 
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Proof. By (12.81) and Ito isometry, it follows that 


miP-pTL^in,)] = E[\\ / Gix,y)dW{y)-Gix,y)dW>^{y)\\i.^^^^] 


= E[ \ G{x,y)dW{y)-G{x,y)dW\y)\^dx] 

J Qfe ^ f 

n N p TT-l p p 

= E[ 1^/ G{x,y)dW{y)-\\K,\\-^Y. G{x,z)dz dW{y))\^dx] 

O-b i—\ i —1 Ki 

K,\\-\G{x, y) - G(x, z))dz dW{y)\^dx\ 

\G{x,y) — G{x,z)'[^ dxdzdy, 


(2.14) 


» N N 

N N ^ 
i=l i=l 


'Ki JKi 


iKi JKj Jnt 

where positive constant C depending on sup, \Ki\. Thus, Lemma yi?f \ leads to the completeness of the 
proof □ 

Theorem 12.31 implies that converges to p as h —?> 0 in a bounded domain Hf,. In Section 3, H;, 

will be choosen as a computational domain. 


3. PML model for convected Helmholtz equation 


Let us study the modified PML method which is described in Section 1. We first prove that the 
solution of PML model converges to the solution of (12.41) in Hf,. Then, it follows from combining with 
the result in Theorem 12.31 that the stochastic error estimation between solutions of (jl.ll) and PML 
model. 

As in [3], the presence of inverse upstream modes should be considered to avoid an unstable 
solution. The presence of instabilities was studied in [S] by using group velocities. We recall the 
modified PML model as follows. The problem (12.31) is reformulated by 


- EH 




)y 


"'dxi ' 1 — 

We now apply the following substitution [2 

d Mki 


d'^ph 




1- A/2 


p'^ = f + W^ 


n. 


+ 


, .,9 Mki 


— lUJ 


dxi 1 — M2 

Here, the complex function a{xi) is given by 

a{xi)= . , / V 

—lUj + a[xi) 

where oj > 0 and cr(xi) is a positive real-valued function in C^(IR). 

Example 3.1. Let 


0 

r(xi) = I CT- 


(3.1) 


(3.2) 


(3.3) 


(3.4) 


if Xi e {x~ ,x+), 
tT+ • (xi — X’*')^ if Xi > X+J 

tT_ • (x“ — Xi)^ if Xi < x~ , 

where a± > 0. This type of u is usually applied for acoustic wave problems. Here, constants a± are 
associated with the reflection coefficient. Refer to [laiiH] for details. 


At this point, we have considered infinite PMLs, {x G xi > x+} and {x G H| xi < x“}. 
In practice, we need to build bounded layers of the finite length L. Thus, let us define two PMLs, 
= {x G x+ < Xi < x+ -I- L} and = {x G x“ — L < xi < x“}. Recall that we defined 
the computational domain Hf, in Section 3. Denote that :={xGH|x“—L< xi < x+ -I- L} is 
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the truncated domain and := {x G ^l\ xi = x^ ± L} is the external boundaries. For simplicity, we 
consider the Dirichlet boundary conditions on Let be the solution which satisfies 


dp^ 

—— = 0 on^nn^sii, 

0 x 2 


5V 

dxl 


e 


krh 


l-IvP 


p^^ = f + W 


in rir 


(3.5) 


= 0 onE±. 

Let Vl = {qG H^{V,l)\ q = 0 on E^}. A variational formula of (13.511 is written by 


where the form aQL{-, •) is 


with 


ani(p,g) = -/ {f + W)qdx, 'iqGVL, 


anr- (p, q) = bciL {p, q) + Cql {p, q) 

1 fc2(l-a2M2) 


a{xi) 1 — 


pq dx, 


bQ.L{P^q)= / [Ma^p) -Vqdx - / 

J J 

CQ,^(j),q) = —2ik / a(a;i)V • ( \ q dx — ikM / a'{xi)pqdx, 

Jql \ ^ j Jcil 


(3.6) 

(3.7) 

(3.8) 


where Ma = 


{1-M^)a 0 

0 i 

a. 

Now we prove the following theorem. 

Theorem 3.1. The problem US. Oil is of Fredholm type. 

Proof. The idea of the proof is based on [4] and [8]. First, we check that there exist a bounded operator 
Kl on such that 

1 kf{l — a^M^) 


{KlP, q)m{Q.L) = cql {p, q)- ( 

Jo,. 


siL «(2;i) 1-M2 


+ l)pqdx, Vp, g S i7^(n^). 


By the compactness of the embedding of into the operator Cl is compact. Let us 

define sn^{p,q) = aQL(j),q) — {KLP,q)H^{p,^). Taking the real part of sni(g, g), it is derived that 


R-e(sni (g, q)) = ( ( (1 - M‘^)Re{a) 

JO.L \ 


dq 


dxi 


+ Re(-) 


dq 


> inf 


(1 - M2)a;2 


dX2 

2 


+ \q? dx 


|Vg| da; +||g||i 2 (n^) 


(3.9) 


>C||g||? 


EGnL \cr2(a;i) + P ^ 

where C = infa;gni This implies that |sni,(g,g)| > C'lkllffi(ni,)' Thus, applying Lax- 

Milgram theorem and Riez representation theorem, there exist a bounded operator Sl on 
such that 

{SLP,q)H^(nc = ^^l{p,(i) 'dp,qG 

and Sl is of Fredholm type. Therefore, from Corollary 4.47 in [I], Sl + Cl is of Fredholm type. □ 

Since (13.61) is of Fredholm type, the problem (13.61) is well posed if and only if the homogeneous 
problem has only the trivial solution p = 0. In order to ensure the well-posedness, it will be discussed 
in Section 3.2. 

Next, we discuss reduced problems which are equivalent to (12.3p and (j3.5ll . These reduced prob¬ 
lems have two main advantages: first, they provide a result of existence and uniqueness of the solution. 
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Second, numerical methods can be used to solve them, since they are posed in a bounded domain fit,. 
In this setting, the error analysis of the PML method is established in a bounded domain 


3.1. Reduction to a bounded domain 

Let us reduce the problem (12.31) in H to a problem in a bounded domain for the error estimation. The 
idea is based on that in [4]. To derive the appropriate boundary conditions for the reduced problem, 
we first consider the homogeneous problem of (12.3p . i.e. the right-hand side is 0. In this case, the 
solution of (j2.3p is represented by 

(3.10) 


where n G {0} UM, 4>o{x2) = 4>n{x) = y | cos(^^^^^), n = 1,2,3, • • •. Here, a;i-axial wave numbers 

satisfy that 

-{I- -2kMp + k^ = (3.11) 

Let Kq — —^==. Let us denote that [iLo] is the integer part of Kq. We set [Kq] = Nq. Then a;i-axial 
wave numbers are given by 


-kM± 


1 - M2 


, if n < No, 


-kM ± M^) -k^ 

^ -, if n > No. 


1 - M2 


(3.12) 


The solutions (x) or (cc) are outgoing waves in the direction of cci —>■ oo {xi —>■ —oo, respectively). 
If P^ is real, then are called propagating modes. The group velocity is positive for p+ and 
negative for p~. The number of propagating modes increases as the Mach number M —>■ 1. In the 
convected Helmholtz equation, it is well-known that inverse upstream modes p^ (x) which have a 
positive group velocity and a negative phase velocity ^ is appeared in [4] . If P^ is complex, then p^ 
is called an evanescent mode. Evanescent modes p^ or p~ are exponentially decaying when xi —>■ oo 
(xi —>■ —oo, respectively). 

To construct the posed problem from (12.311 . some preliminaries are concerned in the following. 
Let us choose a bounded domain Ub = {x € H, x~ < xi < cc+j which contains the support 
of f{x) + W^{x). Define are two boundaries {x G D| cci = x^}. The inner product (•,-)i2(s+) 
(respectively (•,-)l 2(E-)) is defined by 


(M,u)i 2 (s^)= / uv dx 2 . (3.13) 

dE+ 

We now discuss the reduced problem in fib which is equivalent to (12.3p . According to [4], the 
problem (12.311 is reduced to 

+ If ^ ^ 

1^ = 0 ondfindflb, (3.14) 

0X2 

9Pb _ rr, h ^ ^ 

T±Pb ^±) 

where the vector n denotes the outward unit normal to Here, the Dirichlet-to-Neumann(DtN) 
operators T± : i7^/^(E±) —>■ are defined by 


t±W = tJ 2 'iPtii^,(l>n)L^{S±)(l)n{x2). 

n—0 


(3.15) 
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The DtN operators T± play as exact nonreflecting boundary condition so that two problems (12.31) and 
(13.141) are equivalent. Furthermore, a weak from of (13.141) is written as follows: 


andPb^9) = - [ if + W'^)qdx, Vg G 
J Qb 

where the sesquilinear form aa^(-, •) is 

an„ (p, q) = bn^ {p, q) + Cn, {p, q) 


(3.16) 


(3.17) 


with 


bnt{P:q)= / / {MVp)-Vqdx+ / k^pq dx + {T+p,q)Y,++ {T_p,q)Y,_, 

J Qf, J Cl J Cl 


(3.18) 


Here, the brackets (•,-)e+, are the natural duality pairing of and H ^/^(S±), re¬ 

spectively. The well-posedness of the problem ()3.14|) is served in the next Theorem. 


Theorem 3.2. We assume that 

k ^ v'l - M2 ^ (3.19) 

for n G IN. Then the problem i2.4\) is well posed. 

Proof. Refer to Theorem 2.2 in [3]. □ 


We now start with reducing the PML problem (13.5|) . After doing this, the error analysis can be 
derived from the reduced problem. More precisely, our aim is to compare the solution p^ of (13.141) 
with the solution p^ which satisfies: 


Mki 


=0 on n d^lb, 

0X2 

dP^ rrL L V 

-^=-T±Pb onS±, 


d^Pb 

dxl 




1-M2 


Pb = f in Hfc, 


(3.20) 


where the vector n denotes the outward unit normal to . Here, the DtN operators : 77^/^ (E±) — 
i7“^/^(E±) are defined by 


^±(V’) =T'^iVn{'^{xi,-)An{-))L'^(Y.±)(t>n{x2), 


n—0 

= Pi-- 


1 — e 


■^n T'n 
+ _ /Q- ^ 




0 a(x+ + s) 


ds 


Pt-Pn 


1 — e 




0 a(x--s) 


ds 


(3.21) 


Let us prove that the solutions of problems (13.51) and (13.201) are equivalent in Hf,, i.e. = p^ 
Then by using p^ and p^, we can obtain the error estimation with respect to p^ and p^ in Tib. The 
next Theorem asserts that = Pb ■ 


Theorem 3.3. Suppose p^ and p^ are satisfied LS. ,51) and liS.20\) . respectively. Then p^la^^ =Pb- 
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Proof. The key idea is to find a correct boundary condition on In the solution satisfies 
a homogeneous equation. Therefore, we can apply the separation of variable method. Let us consider 
the solution in the domain We define 


Then satisfies 


i’nixi) = e" 


, d Mki 

a{xi){— + 


dxi 1 — 



(3.22) 

Mk 

■)^tixi) =*(/3n + 1 _ ]p[ 2 ')'^n{xi)- 

(3.23) 


It follows from (j3.23|) that 'ip^{xi)(j)n{x 2 ) satisfies the homogeneous equation problem of (j3.5l) with 
boundary condition on dVl fl d^b- By the Dirichlet boundary condition on the solution p^ in 
is given by 


P+{x) = ^(p+(a;+,-),(/>n(-))s+(A+i/’;!:(a:i) {Xl))(j)n{X2), 


n=0 


where 


jPi /c 


An =T- 


0 a(:r++s) 


ds 


e 

rL 1 


fo a(x+ + s)^^ _ J^O a(x++s) 


ds 


(3.24) 


(3.25) 


By (13.191) and (13.3L (/3+ — j3.^ ) a(x++s) ^^ ^ positive imaginary part ds. Therefore, 

it can not be a element of 27rZ := {27r2;| 2 ; € Z}. In other words, the denominator of Afi is not zero. 
In this reason, and the formula (13.241) are well-defined. 

Since Q;(a;i)|i;_^ = 1 and a(xi) is continuous, an exact boundary condition on E_|_ satisfies that 

Niki 


dPb I / d Mki 

= »(-.)(— + 


dxi ' 1-M^ 

Mk 


Pbk+ 


= ^(p+(a;+,-),(/>»(-))s+(*A+(/3+ -k (/?» + j 

n—0 

-f^{p^{x+,-),M-)h+(.iA+ + lAf 


Mk 


- M2 
)(t>n{x2) 


))</>n(a;2) 


(3.26) 


n=0 

00 


= i^{Pb{x""r),^n{-))T.+ {A+l3t + A^(3^)(j)n{x2). 


n—0 

Repeating the argument for the boundary condition on E-|_, an exact boundary condition on E_ is 
dPb 

o—IS- = 

OXi 

n—0 

where 

e 


|s_ = -i^{p^{X-,-),<l)n{-))i:+{B+Pt + BnPn)(t>n{x 2 ), 

fo 


(3.27) 


Bt=T 


2 .{X--S) 


-ds 


ipt If' -rds iPn ff —A-r 

g a(x —s) _ g a(x —s) 


ds 


Let i/+ = A+I3++ and = (H„ /3+ -h H+/3„ ). Then, we complete the proof. 

A weak from of (13.201) is also written by 

ahtiP^<l) = -[ {fo + W^)qdx, Vg G M^(Ht,), 

where the sesquilinear form aQ{-,-) is 

aL {p. q) = bh, (p, q) + Ck (p, q) 


(3.28) 

□ 

(3.29) 


(3.30) 
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with 


q)= [ Vp) -Vqdx- [ P pq dx + {T^p, g)s+ + {T^^p, q)^_, 

^ Q,L 

ch,{p,q) = 


(3.31) 


Remark. Suppose ip^{xi) are outgoing propagating modes, i.e. n < Nq. Then ip^{xi) are modes 
reflected by Amplitudes of the reflection coefficients i?" for each 'ip~{xi) are 


RZ = 


A+ 


’ 1-M^ 




^-ds 


For evanescent modes V'n (a^i) (n- > A^o), amplitudes of reflection coefficients are 


(3.32) 


RZ = 


A + 
■^r). 


= e 


2kL 

1-M'^ 


r-1 


(3.33) 




)ds 


T 2fc Q 

In other words, reflection coefficients of waves caused by Y-'^ are bounded by e ^ , 

where C is a positive constant depending on Kq. For example, if a is that of the example 13.11 then 
upper bounds of reflection coefficients are 


2C jS 

3cq(1-M'^) 


and 




(3.34) 


for propagating modes and evanescent modes, respectively. By using (13.321) and (13.331) . we check that 
the solution p^ converges to the solution as L —^ oo in the following. Suppose the length L of layers 
goes to infinity. Then amplitudes of reflection coefficients converge to 0 by (I3.32|l and (|3.33|1 . This 
means that all modes 'ijj~{xi) reflected by PML are exponentially attenuated. The rigorous proof of 
the convergence will be studied in the next subsection. 

Furthermore, the well-posedness of the problem (I3.20|) follows from Theorem [3T] and (13.2211 . The 
well-posedness of the problem is shown in the next theorem. 


Theorem 3.4. Suppose that Ii3.19\) holds. Then the problem h3.2(A) is well posed. 

Proof. Applying Theorem 13.11 it is enough to show that the homogeneous problem of (|3.20|1 has 
only trivial solution. By applying (I3.22p in Theorem 13.31 and replacing x'^ with 0, the solution of the 
homogeneous problem is decomposed as follows: 

OO 

Pbi.x) = + D~if~{xi))(t>n{,X2), (3.35) 

n—0 

where are complex constants. The boundary condition = —T^p^ on Y^ implies that 

Dt{v^-Pt)=D-{vt-Pu) = ^ forn = 0,l,2,.-.. (3.36) 

Since — Pt) = ~ Pn) = 0 equivalent to (/3+ — Pn) =0, the assumption (13.191) leads to 

=0, for n = 0,1, 2, • • •. Consequently, p^{x) =0. □ 


3.2. Convergence and error estimates in flf, 

Recall that the object of our study is to establish the error analysis with respect to p of and p^ 
of (j3.5p . As discussed in section 2, we applied the discretized white noise for building more regular 
problem (j2.3|) and provided the error estimation E[\p — Furthermore, it was shown from 

()3.14p and Theorem 13.31 that p^ and p^ are equivalent to p^ and p^ in fib, respectively. Therefore, we 
would be able to obtain the error between p^ and p^ which leads to the final result in Theorem 13.71 In 
this reason, we first study the error analysis with respect to and . To do this, some preliminary 
observations are discussed in the followings. 
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Let us consider and which satisfy variational formulas: 

ifo + W^)qdx, VqGH^iflb), 


(3.37) 


o-htiPb ^ ^) = - [ {fo + W’^)qdx, yqGH^iQb), 

J Qfo 

respectively. We also define a function d{xi) = (1 A In this setting, the next lemma is proved. 

Lemma 3.5. Suppose the length L of layers is large enough. Then there exists positive constants Ci 
and C 2 such that for all p,q € H^{flb), we have 

Precisely, C 2 = t^( 1 A _ f). 

Proof From (12.41) and (13.6L it follows that 

la-ndP^q) - an,ip,q)\ < \{T+p,q)^^ - {T^p,q)^+\ + \{T.p,q)^_ - {Tfp,q)^_\. (3.39) 

Since the estimation of two terms in the right-hand side of (13.3911 are analogous, we only derive the 
upper bound of the first term \{T+p,q)-s^. — (P+P, (z)e+|- Let ( = p\s+ and 77 = gls.,.. From (13.151) and 
(|3.21|) . it is derived that for C G i7^/^(S+), 


(r+ - )C = - ^ 7(/3+ - l^+KnM^2), 

n—0 

where Cn(a^i) = (C(2^ii'): </^n(-))L^(s+)- Thus, it is shown that 

00 

((T+-Tf)C,??)E+ = 

n—0 

where rjn = {p,(fn)L‘^{Y,+)- Hence, (13.3p and (13.211) lead to the inequality 


(3.40) 


(3.41) 


K(r+ - T^)C,r;)Ej < \Pt - K\\Cnnn\ = E 


- Pn 


n—0 


^ h _ eh/3n-/3+)(L+*/<f ^d£±±£)d.) 


'ICra^ral- (3.42) 


On the other hand, if z G C, Im(z) < 0, and |Im(z)| is large enough, then the following inequality 


holds: 


|1 _ e“| > - 1| > 


To apply the inequality (13.431) for the term |1 — 
cases in the following: 


L <T(x++a) 


(3.43) 

in (I3.42F we consider two 


I. propagating modes : suppose /3+ — /3„ = y^l — i.e. n < Nq. Then, Im((/3„ — Pn){L + 

i ^1 — ^ Jq Clearly, this quantity is negative. By choosing large 

enough L, the inequalily (13.431) holds for propagating modes. Therefore, it holds that 




4/c 


^2 2fc 

^ 1-M^ 


< 


1-M2 

4fc 


' 1 fo 


1 P 


-ds) 


2 ^ l-L fl£l±£)^g) 


(3.44) 


1-M2 

H. evanescent modes : suppose — Pif = 




1 i.e. n > A^o + 1- Then, Im((/3„ - /3+)(L -|- 
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i ^^ds)) = — \J^ ~ choosing large enough L, the inequality (13.431) is verified for 

evanescent modes. Thus, 


\fit-K\ < — 


4fc 


M2y Kl 


< 


From (j3.44l) and (I3.45p . we obtain 


n=0 


Nq 


l-M"^ Kn 


2k rL cr(a:++5) 


— 1 e 

1-M^ Y 

2kL / 

(Nn + i)2 , 

Y 

0 

‘^"^\CnVn\ 


4fc 


E 


^ 1 — m'2 


2kL / (^0 + 1)^ 


1 - M2 ^ Ko 

n=iVo + l 


l-M^ V K, 


iCnVn 


n =0 ^0 


<C-e 

-1 - M2 


(3.45) 


(3.46) 


According to the trace Theorem in [16], it follows that 

|((T+-T^)p,q)sJ<Cie-^^/o"-("^+*)‘^1b||^qn.)lkll/rHn.). (3-47) 

By repeating the procedure for (I3.47F a similar result for |((T_ — T^)p,q)s_ \ is shown as follows: 

|((r_ -rf)p,<7)E_| < (3-48) 

□ 


Now, we turn to the error estimation \\p^ — ^ ~ i7^(f2b). Then linear operators 

A and in C{V,V') are defined by 

{Ap,q)v' = an^{p,q), 

r r 3.49 

{A^p,q)v' = aa^{p,q), 

for all p,q G V. Here, the operator A is the same as that of p.424 in |4]. By Riesz representation 
Theorem, A and A^ are uniquely determined. Moreover, operator norms ||.4||£(v_y/) and \\A^\\c{v,V') 
are bounded by upper bounds of forms anfc(-, •) and aQ^{-, •), respectively. It follows from Lemma [37^ 
that 


M-^^ll£(y,v') < + (3.50) 

Since p^ and p^ satisfies that for all g S F, 

{APbyq)v'=-{Lq)L^{n), 

. (3.51) 

(-4 Pi ,g)y'=-(/, g)L2(n), 

respectively, it holds that 

(-4^(Pb - Pb), q)v' = {{A^ - A)Pb, q)'v (3.52) 

for all qGV. 

From (I3.50p and (|3.52|) . the error estimation for |bb — P^Wh^^o.^) is established as follows. 
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Theorem 3.6. There exists a large positive constant L such that for all L > L, is an isomorphism 
on H^{flb). Moreover, we have 

7’+J 


(3.53) 


where C\ and C 2 are positive constants depending on k and M. 

Proof. The argument of the proof is similar to that of Theorem 4.4 in [4]. The operator can be 
rewritten as A + {A^ — A). According to [1] , A is an isomorphism on the Hilbert space V . Therefore, 

{Ap,q)v> = {p,A-^q)v. (3.54) 

We consider the problem: finding u €V which satishes 

{A^u,q)v' = {g,q)v', 

for g gV and all q gV. From (13.541) . the problem (13.551) becomes 

((/ + A~^{A^ - A))u, q)v = {A~^g, q)v- 

From (j3.50L it holds that for L > L, 

This leads to 

\\A-\A^-A)\\<1, 

where || • || is a formal operator norm. By the Banach fixed point theorem in 
admits a unique solution. Furthermore, it follows that 

1 


(3.55) 

(3.56) 

(3.57) 

(3.58) 

the lineear map A^ 


\\iI + A (A A)) II < 

By (13.561) and (13.591) . the following inequality holds: 

II II . M~^g|lv 

l-\\A-^{A^-A)W 

Let us set u = p^ — p^. Then, by (I3.52|) and (I3.60L we obtain 

\\A-HA^-A)pnv 


\\Pb Pb IlffHn^) < _ ||_4-i(^L _ _4)|| 


< 


\\A-mA^-A\M\\v 


(3.59) 


(3.60) 


(3.61) 


1-\\A-^{A^-A)\\ • 

Since \\A^ — A\\ is small enough by choosing some L, it is possible to taking a positive constant C > 1 
as a upper bound of ■ Consequently, the proof is completed from p.50l) . □ 

Theorem 13.61 implies that p^ converge to in Hf, as L —>• 00 . Applying two estimations (|2.13|) 
and (I3.53p . it follows that the next error estimation holds. 


Theorem 3.7. Suppose that S 8.191) holds. Let the length L of the layers is large enough. Then there 
exist positive constants C\ and C 2 depending on k, sup \Ki\ and M such that for 0 < e < 1, we obtain 

E[\\P - P^Wh^n,)] < Ci{h^-^ + ^0 A=++s)ds ^ ff a{x--s)dsy 

Proof. By the triangle inequality, it follows that 


E[\\P - P^WhiaJ =E[\Kp - /) - (/ - P^)\\Un,)] 

<‘^E[\\p - p^||i2(nj)] + 2i?[||Pb “ Pb\\L^(nh)] (3.63) 

Therefore, by Theorem 12.31 and Theorem 13.61 the next inequality is proved: 


(3.64) 
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□ 

Theorem O asserts that the solution of the problem (13.611 converges to p in Ilf, as L —?> oo 
and h —>■ 0. Therefore, we obtain an approximated solution of p by solving the variational problem 

(1^ . 


Conclusion 

We have studied the stochastic convected Helmholz equation in an infinite duct. Since the regularity 
of the solution of SPDEs is generally weak, an alternative problem has been constructed by using the 
discretized white noise. In this setting, we have proposed the modified PML model which omits the 
presence of inverse upstream modes. Applying modal expansion approach, an error analysis of the 
PML model has been provided. Finally, the stochastic PML-truncation error of the solution has been 
established in the sense of • jHa]- 


Appendix 


The proof of Lemma 12.21 

Proof. By (j2.6ll and (12.7L we obtain 

OO 

\G{x,y) - G{x,z)\ <2|$(x, 2/) - $(a::, 2 :)| + ^ \<^{x,et, + y) - $(a;,e+ + z)\ 

n—1 

OO 

+ XI ^n+y)- ef + z)\ 


n—1 

OO 




X, e+ - z) 


(3.65) 


n—1 


n—1 

=: 2/o + Ii + I 2 + Id. + li- 


X, e„ - z 


To prove the inequality (I2.12L we will derive the following inequalities: 

I if dx < C\y — z\^~'^, f = 0,1,2,3, 


(3.66) 


where C is a positive constant depending on M and k. 

First, the case of /q is considered. From Lemma o and the boundedness of the domain Hf,, it 
follows that 


/ lfdx<C{ \\n p{x — y) — \n p{x — z)\^ dx + / \V{x,y) — V{x, z)\^ dx) 

^G j \\n p{x — y) — \n p{x — z)f dxC\y — zf. 

J Qb 

For 0 < e < 1, the integral | lnp(a: — y) — \np{x — z)\'^ dx is reformulated by 

[ I \np{x -y)- \np{x - z)Y{p{x - y) - p{x - z)f-^ ( [ 

Jn,, \Jo 


(3.67) 


ep{x - 2 /) + (1 - 9)p{x - z) 


2-e 

d0 1 dx. 


(3.68) 
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By the triangular inequality and Holder inequality, (13.681) is bounded by 

|lnp(a:- 2 /)-lnp(a;-z)|" ( / — -— ^ - -do) dx 

Jn, \Jo dp[x-y) + [l-d)p{x-z) J 

<C|y-zp"'^/' |lnp(a;-y)-lnp(a;- 2 ;)|" 

Jn, \p{x-y) p[x-z)J 

/ c \ / /* 1 1 \ 

<C\y-z\‘^-''\^J^\\np{x-y)-liip{x-z)\^dxj I [a- _ ^| + [a- _ ^| l^^dxj 

(3.69) 


Since p{x — y) and p(x — z) are bounded by some C > 0 for x, y, z S Hf,, the integral | In p(x — y) — 
In p{x — z)\^dx has an upper bound 

rC 


c f 

1 lnp(x 

J Q^b 


f 1 1 

1 1 1 

JQb 1 \x-y\ 

c [ 

1 1 

\x-y\ 

J fib 


3(2-e) 

3-' dx is also bounded by 


3-' dx < C / r 3-e dr < oo. 
Jo 


(3.70) 


(3.71) 


2-e 


(3.72) 


By (|3.69|) . (13.701) . and (13.711) . we conclude that 

[ lQdx<C\y-z\ 

J Oj, 

where C is a positive constant depending on M and k. 

On the other hand, Hankel function Hq(z) has an asymptotic behavior 

iji(z) = /4)(i ^ 0(1)) (3.73) 

V TTZ Z 

for large real z. To prove if dx < C\y — zp“'^, we will use (|3.73|) in the following. We choose a 
large uq > 0 such that p(x — y — e^) and p{x — z — e+) have a lower bound C > 0. Then, $(x, e+ + y) 
and $(x, el + z) have an asymptotic behavior p.73l) for all n > uq. Choosing sufficiently large no, it 
follows that 

\^{x,e+ +y)-$(x,el +z)| < \H^{kp{x - y - e+)) - H^{kp{x - z - e+))\ 

<C\- ^ ^ ' 


\l P{x-y- el) p{x- z- el) 


< C 

< c 

< c 


\p{x -y-e+)- p{x- z-e+)\ 


p{x - y - el)p(x - z - et){\lp{x - y - el) + \/p{x - y - el)) 

\p{y- z)\ 

rfid^ 

\y-z\ 

rfid^ ’ 


(3.74) 


where C is a positive constant depending on M and k. The inequality (|3.74p implies that 

OO OO ^ 

|$(x,el+ y) - $(x,el+z)| < Y - ^1- 


(3.75) 
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For n < no, the same procedure as in the case of Iq is available. Hence, we obtain 

/ \^{x,e^ + y) — + z)\'^ dx < C\y — for all n < no- 


By (13.751) and (13.761) . for small h > 0, it holds that 

rto —1 


P /^u J. « 

/ lldx<Ci^ / |$(a;,e++2/) - $(x,e+ 

Hfe n—'V 

P oo 

+ C2 I ^ ($(x,e++ ?/) - $(x,e++ z))|^dx 

n—no 

< Ci\y - z\'^~^ + C2 [ \y-z\'^dx 

r2>i 


(3.76) 


(3.77) 


< C\y - ", 

where y,z € Ki and C is a positive constant depending on k, d, |Hb|, and M. 12 , 3,4 cases are proved 
by the same method as in the case of Ii. 

□ 
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